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We give a heuristic argument for disorder rounding of a first order quantum phase transition into 
a continuous phase transition. From both weak and strong disorder analysis of the the A^-color 
quantum Ashkin- Teller model in one spatial dimension, we find that for A'^ > 3, the first order 
transition is rounded to a continuous transition and the physical picture is the same as the random 
transverse field Ising model for a limited parameter regime. The results are strikingly different from 
the corresponding classical problem in two dimensions where the fate of the renormalization group 
flows is a fixed point corresponding to A'^-decoupled pure Ising models. 

PACS numbers: 



The effect of disorder on continuous classical phase 
transitions has been intensely studied over many 
decades, [T] but less is known about its effect on first order 
transitions. Imry and Wortis 2] argued that arbitrarily 
weak disorder can actually round a classical first order 
transition into a continuous transition. Subsequently, 
Hui and Berker [11 and Aizenmann and Wehr 0|, Q 
have made important contributions to this topic. In 
light of numerous recent experiments involving quan- 
tum phase transitions in systems that inevitably contain 
many sources of disorder, the issue of disorder round- 
ing in first order quantum phase transition (QPTi) has 
acquired considerable urgency. 

There are two important questions: (a) Can disorder 
convert a QPTi to a continuous one? (b) Might the theo- 
rem for a disordered classical system require modification 
in a quantum context? There is, as we shall see, a simple 
intuitive affirmative answer to (a), but (b) is more subtle 
because statics and dynamics are entangled in a quantum 
phase transition. To answer (b) and to explore more fully 
the issues involved, we consider below a model and study 
it in considerable detail. It is important to note that for 
even the classical random bond Potts model for q > 4 
and d = 2, for which the pure system has a first order 
transition, the critical exponents of the disordered sys- 
tem do not belong to the simple universality class of the 
pure Ising model — although the critical exponent « 1 , 
all other critical behavior is different Q. 

We will answer question (a) by a heuristic argument 
patterned along an argument by Berker Q , although we 
differ in our analysis for the case of continuous symme- 
try. Let the disorder couple to the Hamiltonian in such a 
way that its symmetry is unchanged. For example, dis- 
order may couple to a nearest neighbor bond (generally 
to energy-like variables) without affecting the symmetry. 
In contrast, a site random field breaks the symmetry ex- 
plicitly. Consider changing a generic tuning parameter 
g that may be the ratio of the amplitudes of two non- 
commuting terms in the Hamiltonian, which controls the 
quantum fiuctuations and results in a first order transi- 



tion. 

We would like to show that in the presence of disorder 
coexistence of phases is not possible at the transition, 
and the quantum fluctuations do not have a scale. If 
this is true, and if the state corresponding to g = is 
still a broken symmetry state (this is why we imposed 
the specific requirement on disorder earlier) and the g — 
oo is the quantum disordered state, the conclusion must 
be that the transition has been converted to a quantum 
critical point. 

The proof is by contradiction. Assume that the first 
order transition is at gc, implying coexistence of phases. 
However, in the presence of disorder there will be lo- 
cal fluctuations of gc- Thus, within a putative quantum 
disordered region, randomness can nucleate an ordered 
region of linear dimension L, with a gain in the volume 
energy (x (assuming central limit theorem), while 

the price in the surface energy is oc L'^~^. The same is 
true for a putative ordered region. Therefore, for d < 2 
(discrete) , the picture is that of a "domain within do- 
main", and there is no scale, as required for a QPTi. In 
contrast, for g < gc nucleation of one broken symmetry 
phase within another does not gain any energy (disorder 
does not break the relevant symmetry), but the surface 
energy is increased. Therefore, by contradiction, coexis- 
tence of phases is not possible, and the transition at gc 
must be continuous. We have tacitly assumed that the 
transition involves a broken symmetry. If this is not the 
case, there is no particular reason for a sharp transition 
to remain at gc, and the disorder will simply smear out 
the transition. 

The case of continuous symmetry is a little subtle. 
While one may like to argue that the domain wall en- 
ergy is L"*"^ [7|, as in the Imry-Ma argument, this is 
typically incorrect. If the domain wall connects a broken 
symmetry state with an unbroken symmetry state, where 
the amplitude of the order parameter vanishes, the do- 
main wall energy is still L'^~^, as in the case of discrete 
symmetry. However, if the QPTi is effected by tuning a 
"magnetic field" that changes the state from one broken 
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symmetry direction to another, the domain wall energy is 
indeed L'^~^, and the borderline dimensionality is d = 4. 
Of course from Mermin- Wagner theorem there is no long 
range order in d = 2 at any finite temperature, by either 
first order or continuous transitions; so, for the classical 
case the question is moot at d = 2. 

There are no simple arguments known to us for the bor- 
derline dimensionalities, but from the rigorous version of 
the Imry-Ma argument for the random field case 0, H, S] , 
it is safe to conjecture that the above argument should 
also hold for these cases, because there is a close con- 
nection between the two problems as noted by Imry 0. 
Note that the dynamic critical exponent z does not enter 
in this analysis — all we need is the extensivity of the 
ground state energy and its normal fluctuations in the 
thermodynamic limit. The principal disordering agent 
that washes out the coexistence is the fluctuations due 
to impurities and not quantum fluctuations. Quantum 
fluctuations can only help the process of smoothing the 
coexistence. Of course the fate of the system in dimen- 
sions higher than 2 must depend of the quantum fluctu- 
ations. In addition, the actual dynamics of the system 
must involve these fluctuations as well. 

It is useful to note that, in contrast, the Harris crite- 
rion [l[ that determines the influence of impurities at a 
critical point, inter alia a quantum critical point, does 
depend on the quantum dynamics. In order to see this, 
let us rephrase the Harris criterion along an argument es- 
sentially due to Mott [10,] . On one hand, disorder in a do- 
main of linear dimension, ^, the correlation length of the 
pure system deflned by the decay rate of the equal time 
correlation function, will give rise to fluctuations of the 
quantum critical point gc of fractional width Ag ~ C 
On the other hand, Ag must be less than the reduced dis- 
tance from the quantum criticality implied by ^, that is 
^ ^-i/iy^ ^jj^g criticality to remain unchanged, where 
v is the correlation length exponent of the d-dimensional 
quantum system at zero temperature. Hence, > 2/d. 
Otherwise, the system may be described by a new disor- 
der flxed point for which the same relation will apply with 
the replacement of the critical exponent of the pure sys- 
tem by the critical exponent of the new fixed point, as in 
the theorem of Chayes et al . In either case the quan- 
tum dynamics is important because the relevant length 
scale close to the critical point is the diverging correlation 
lengthy ^. By contrast, the argument involving QPTi is 
restricted by a finite correlation length, hence the balance 
is between the volume energy and the surface energy of 
a fluctuating domain nucleated by the impurities. 

We now turn to the question (b) and study, using both 
a perturbative renormalization group and a real space 
decimation procedure [13, [13, [3|, ^ one-dimensional 
quantum s pin chain, the random A'^-color Ashkin- Teller 



model 15, 3, [17| in the presence of disorder. We con- 
sider the regime in which the pure model has a QPTi. 
The corresponding classical problem in two spatial di- 



mensions where the quenched disorder is isotropic, the 
renormalization group flows curl back to the pure decou- 
pled Ising fixed point 18, 13], at least for weak couphngs. 



The Hamiltonian of the system is [2C 
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Here, Latin letters index lattice sites and Greek letters 
label colors and the cr's are the Pauli opertaors. The Ji 
and hi are random variables taken from a distribution 
restricted to only positive values, while e is a disorder 
independent positive constant. For the random trans- 
verse field Ising model (RTFIM, e = 0), a local gauge 
transformation may be performed to make all couplings 
positive, so the original couplings can take negative val- 
ues. The coupling between colors destroys this freedom 
so we must restrict ourselves to positive couplings from 
the outset. We have parametrized the system so that 
the intercolor couplings are proportional to the bond or 
field at that site. We also restrict ourselves to e > 0. 
Note the invariance of the Hamiltonian with respect to 



the following duality transformations: a\ ,- V. 



J") 
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lJ■l^, 

Ji ~* hi, where /x's are the dual Pauli 
operators. For the uniform system and with N > 3 and 
for e > 0, there is a first order transition from a param- 
agnetic to an ordered state 13, 13, 12, 21). 



For weak disorder and weak four spin coupling, we 
can consider the continuum action in terms of Majorana 
fermions. The random A^-color quantum Ashkin- Teller 
model can be described by an 0{N) Gross- Neveu model 
(GNM) with random mass. [13] The GNM action is given 

by 
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where vp is the Fermi velocity, ^ = J dxdr, and 
^(a) ^ J^(")]t^ = [il;(°')ficr2. In the above equa- 

tion the random mass m{x,T) follows Gaussian white 
noise distribution such that the correlation is given by 
m{x,T)m{x' ,t') = Ai5(a; — x ). After averaging over dis- 
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order using n-replicas, we obtain 
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Here the index a corresponds to replicas. Simple power 
counting shows that A is a relevant operator with scaling 
dimension 1 and g \s a. marginal operator. For even num- 
ber of colors (TV = 2M), this action can be expressed in 
terms of M Dirac fermions and can be bosonized [23| ; for 
TV = 2M + 1 there will be a leftover Majorana fermion, 
which makes the analysis more complex but physical 
answers are the same. A perturbative renormalization 
group calculation following Giamarchi and Schulz[24] up 
to 0(A), 0(Ag) and 0(g2) lead to the following recursion 
relations, where A — Aa, a being the lattice spacing: 
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where z is the dynamic exponent. These recursion rela- 
tions are valid for g/n <C 1. For vf to remain fixed, z 
must vary continuously: 



z = 1 



A , a. 
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There is one unstable fixed point A* = g* = (free 
fermion fixed point with z = 1) and both disorder and 
four spin coupling constants are relevant perturbations 
at this fixed point and flow away to the strong coupling 
regime. This makes it necessary to attack the problem 
using strong disorder renormalization group approach. 

The above flow equations should be contrasted with 
the flow equations of the classical A^-color Ashkin- Teller 
model with quenched disorder. In the classical case the 
the disorder averaged action is local and disorder as well 
as g are marginal perturbations. For this reason one has 
to calculate up to ©(A^), 0(52) and 0{Ag) [H, [H, [H . 
In the replica limit, n — > 0, the flows curl around and 
end up at the decoupled Ising fixed point, at least at 
weak couplings. 

To solve the random system, we will now employ the 
strong disorder renormalization group technique [isl . 
14l . |26[ . There are some similarities with the analysis of 



the quantum gr-state Potts chain for 5 > 4 for which it 
was argued that this transition is described by the in- 
finite disorder fixed point, at least when the disorder is 
strong [2^ . The decimation equations can be inferred 
from the structure of the energy levels. We will demon- 
strate the calculation for a site decimation and note that 
the bond decimation equations follow from duality. If 
the magnetic field on site i is the largest coupling, the 
unperturbed Hamiltonian is 
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and the ground state is 



with energy = 



—Nhi — (2)6/1^. The Hilbert space is written as a tensor 
product of the spins at site i for each color. There are 

N first excited state of the form |— >< > ■ • ■ ), where one 

color has its spin flipped. Each of these states has energy 
Ei^ Eo + 2hi + 2{N - l)h.te. In general, the r*'* energy 
level has (^) states with r colors flipped and energy Er = 
Eq + 2rhi + 2r(N — r)hi€. The coupling of site i to the 
rest of the system is the perturbation part 
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Within second order degenerate perturbation theory, the 
Ising terms connect the ground state to the first excited 
states only, while the 4-spin couplings connect to the sec- 
ond excited states. Thus, the 4^-fold degeneracy of the 
ground state (due to the neighboring spins) is split by V: 



Eo 
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The cross terms yield an effective Ising coupling between 
sites i — 1 and i + I given by 



J 



J^-lJ^ 



h,il + e{N-l)) 
and an effective four spin coupling 

Ji-lJi 



eJ = 



2 h,{l + e{N ~2)) 



(11) 
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Since hi was the largest coupling in the system, the new 
effective two- and four-spin couplings are both smaller 
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than their original respective counterparts. We see that 
the decimation results in new effective bonds and fields 
given by 



J, = 



JiJ i- 



hi = 



hihi 



(13) 



where we have introduced k 
renormalizes as 



1 + {N — 1) e, and e also 



e^{l + {N-l)e) 
2{1 + {N -2)e)' 



(14) 



Equations ^3]) exhibit the duality present in Eq. ([T]) upon 
interchange of the couplings h ^ J. As long as e is 
initially less than some edN), it is clear that e will be 
reduced by this decimation. 

Let n to be the largest energy scale still active in 
the problem and define C = hiil/J, (3 = Infi/Zi, and 
r — Infli/n where ili is the strength of the original 
strongest bond. Then the distributions of logarithmic 
bonds and fields at energy F, denoted by P{C) and i?(/3), 
respectively, satisfy a set of fiow equations. (We suppress 
the dependence of P and R on F.) For P{() we get 



dP 



dP 



+ i?(0) j <P(C')i^(C-C'-ln«) 

+ (P(0)-i?(0))P(C) (15) 



with a similar equation for fiow of R upon the replace- 
ment P ^ i?, as expected by duality. The effect of the 
coupling between colors is to simply to shift the convo- 
lution in the flow equation. If the initial distributions 
are equal, the fixed point coupling distributions are of 
the same form as in the Ising case, and the rescaling of 
energies by F makes the Ins term irrelevant once one 
is below an energy scale on the order of ^i/k [13, 13 ■ 
This is clearly true if e < ec{N) initially. Then e is low- 
ered towards zero, hence k is pushed down towards 1 and 
InK 0. We see, therefore, that for any finite value of 
N, the criticality of the system is that of the infinite ran- 
domness Ising fixed point. Note that from symmetry k 
does not renormalize for either the clock model or the 



Potts model [27[, unlike the present case. So, the above 
analysis is all that we need to perform. 

If the initial value of e is larger than ec{N), n grows 
without bound, so the energy cutoff below which one 
must be to observe universal scaling behavior of RT- 
FIM is driven to zero; in other words, the strong disor- 
der renormalization group analysis breaks down. This 
breakdown of scaling may imply persistence of QPTi 
and will in turn imply an important modification of the 
Aizenman-Wehr theorem. The conjecture is currently 
being checked in numerical simulations f29l| . 

Consider a Hamiltonian, H = Hq + gHi, where Hq 
and Hi commute. A level crossing can take place at 
gc, where an excited state drops below the ground state 



at gc- This will correspond to a first order transition 
and is possible even in a finite system. An example is a 
metamagnetic transition tuned by an external magnetic 
field (soj l . Our work cannot be relevant to this problem, as 
the thermodynamic limit was essential for the argument 
given above regarding the rounding of QPTi by disorder 
into a continuous phase transition. It is easy to see that in 
this case the disorder will merely broaden the transition. 
By contrast, the problem we considered involved non- 
commuting Hq and Hi, and the QPTi was driven by 
quantum fluctuations. 

A striking but simple extension that may also find ap- 
plications to numerous complex strongly correlated sys- 
tems such as organics, heavy fermions, and high- Tc su- 
perconductors is when the QPTi in the pure problem is 
between two ordered states, which from Landau theory 
is generically a first order transition. The heuristic ar- 
gument goes through straightforwardly if we are mindful 
that both sides of the transition involves broken symme- 
tries, albeit of different types. 
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